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ABSTRACT 

Results of numerical simulations of the impact of a common envelope on the matter 
flow pattern near the outflowing component in a semidetached binary system are pre- 
sented. Three-dimensional modeling of the matter transfer gas dynamics in a low-mass 
X-ray binary X1822-371 enable investigation of the structure of flows in the vicinity of 
the inner Lagrange point L\. Taking into account the common envelope of the system 
substantially changes the flow pattern near the Roche surface of the outflowing com- 
ponent. In a stationary regime, accretion of common envelope gas is observed over a 
significant fraction of the donor star's surface, which inhibits the flow of gas along the 
Roche surface to L\. The change in the flow pattern is particularly significant near L 1; 
where the stream of common envelope gas strips matter off the stellar surface. This, 
in turn, significantly increases (by an order of magnitude) the gas flow from the donor 
surface in comparison with the estimates of standard models. 
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1. Introduction 

The observational manifestations of interaction in 
semidetached binary systems (cataclysmic binaries, 
low-mass X-ray binaries, and supersoft X-ray sources) 
are extremely interesting. Comparisons of compu- 
tational results from numerical simulations of mat- 
ter flows in these systems with observational data 
provide a basis for studying and understanding the 
physical processes occurring in these binaries. As a 
rule, numerical simulations of semidetached systems 
have been performed using the generally well-justified 
assumptions that the orbits of the system compo- 
nents are circular and the components' rotation is 
synchronous with their orbital motion. In this case, 
the standard scenario for mass transfer in the system 
is as follows. In the course of its evolution, the donor 
star fills its critical surface (which, under the outlined 
assumptions, coincides with the Roche surface in the 
restricted three-body problem), and mass begins to 
flow through the vicinity of the inner Lagrange point 
Li, where the pressure gradient is not balanced by 
the gravitation force. 

The flow pattern in semidetached systems has been 
investigated using both analytical [1-3] and numeri- 
cal [4-15] models. However, in all these studies, the 
impact of the common envelope of the system on the 
structure of gaseous flows was either not taken into 
account at all, or was taken into account not entirely 
correctly. The morphology of flows in binary systems 
allowing for the presence of a common envelope was 
first considered by us in [16] (henceforth, Paper I). 
There, as well as in this paper, the "common enve- 
lope" of the system refers to the gas filling the space 
between the two components of the system and not in- 
volved in the accretion process (i.e., not belonging to 
the accretion disc). Our results of three-dimensional 
numerical simulations of a low mass X-ray binary 
X1822-371 in Paper I showed the contribution of the 
common envelope to the formation of flow structures 
in the system to be significant. 

The presence of a common envelope also influences 
the flow pattern in the vicinity of the donor star, 
which, in turn, is reflected in the mass exchange rate 
in the system. This problem is extremely important, 
because though the observational manifestations of 
interaction primarily depend on the general flow pat- 
tern studied in Paper I, the evolution of the system is 
determined by the mass exchange parameters. Here, 
we present the results of an analysis based on the 



calculations described in Paper I, but with the main 
accent on a detailed consideration of the influence of 
the common envelope on the structure of gas flows 
near the Roche surface of the outflowing component. 

2. Flow Structure in the Vicinity of L\. 
Results of Standard Models 

Analytical studies of the mass transfer process in 
semidetached systems make it possible to draw con- 
clusions concerning the parameters of flows of matter 
(streams) in such systems. The first detailed anal- 
ysis of matter flows in the vicinity of L\ was car- 
ried out by Lubow and Shu [1], who estimated the 
basic characteristics of the stream using a perturba- 
tion method. In later papers, the stream parameters 
were refined using Bernoulli integral analysis, and the 
dependence of the mass transfer rate on the degree 
of overflow of the Roche lobe by the donor star was 
found [2,3]. The basic characteristics of the streams 
obtained in various studies (including numerical cal- 
culations, e.g., [17]) differ only in details, and are 
currently widely used as standard results in analy- 
ses of mass transfer in semidetached binary systems 
(see, e.g., the reviews [18, 19]). Unfortunately, none 
of these studies took into account the impact of the 
common envelope of the system on the flow pattern 
near the donor star. The aim of our study is to con- 
sider the contribution of the common envelope and 
to obtain new stream parameter estimates using the 
results of three-dimensional numerical simulations of 
the flow in a low-mass X-ray binary system (these 
will be presented in Section 3). However, when an- 
alyzing changes in the flow pattern associated with 
the presence of a common envelope, it is expedient 
to first consider the generally accepted values for flow 
parameters in semidetached systems. 

We will first define the reference frame for our anal- 
ysis. We place the origin of the coordinate frame at 
the donor star center of mass. In our Cartesian co- 
ordinate system (x,y,z), we will direct the X axis 
along the line connecting the centers of the stars, the 
Z axis along the rotation axis, and the Y axis so that 
the resulting system is right-handed. We will set the 
characteristic scale of the coordinate system by plac- 
ing the secondary star at the point (A, 0,0), where 
A is the distance between the two components. This 
means that the center of mass of the system is at the 
point (fiA, 0,0), where = M 2 /(M 1 +M 2 ) is the ratio 
of the accretor's mass to the total mass of the system. 
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It is assumed in standard analytical models that 
the matter at the surface of the donor star (which co- 
incides with the Roche lobe) is characterized by some 
temperature (or local sound speed Co) and density. 
The kinetic energy of this matter (~ c 2 ,) determines 
the degree of overflow of the Roche lobe. The local 
sound speed of the gas for typical semidetached bi- 
naries is much less than the orbital velocity of the 
system components v or it ~ S1A, and their ratio is, as 
a rule, used as a small parameter in analyses. Pro- 
ceeding from consideration of the energetic conditions 
of the gas at the Roche surface, it is not difficult to 
estimate the characteristic dimensions of the matter 
stream that forms in the vicinity of L\. The store of 
kinetic energy of the gas allows it to move in the YZ 
plane passing through L\. Equating the potential en- 
ergy difference and the specific kinetic energy, we can 
derive an equation for the shape of the stream in the 
vicinity of L\, similar to the corresponding equation 
from [1]. Following the adopted scheme for deter- 
mining the YZ cross section of the stream passing 
through Li, we can write: 
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where $ is the force field potential in a reference frame 
rotating with the angular velocity f2 of the system's 
rotation about the common center of mass, in the 
absence of other external forces, excluding gravity. 
The potential in the Roche approximation has the 
form: 
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where the term ((a; — [iA) 2 + y 2 ) is associated 

with the centrifugal force. 

Expanding expression (1) into a Taylor series in y 
and z in the vicinity of L\ and using the condition 

it is not difficult to show that the stream cross section 
is an ellipse with semiaxes 
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where g y (p) ~ 1 and g z {n) ~ 1 arc the coefficients 
of the expansion of the Roche potential in the Taylor 
series in y and z in the vicinity of L\, determined by 
the relation: 

$ = *lx + gyW'y 2 + 9z(p)n 2 z 2 - 

Apart from considering the stream dimensions at 
the point L\, it is also of interest to determine its 
deflection as a result of the Coriolis force. Accord- 
ing to the analysis of Lubov and Shu [1], the stream 
deflection angle is 
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The rotation angle of the stream depends on the mass 
ratio of the components, and, in the simplest case 
(without taking into account the stream's interaction 
with the common envelope gas), it lies in the range 
from -28.37° to -19.55°. 

In more thorough analytical models (see, e.g. [2,3]), 
the Bernoulli integral for the matter leaving the donor 
star surface is used to determine the stream parame- 
ters. The conservation of this integral along a stream 
line, 



u 

T 



7-1 ~ Ll ' 2 ' 7- 1' 

and the assumption that the gas flow velocity is equal 
to the sound speed (the Mach number is 1) can be 
used to derive both the stream dimensions and the 
cross-sectional distribution of parameters across the 
stream. The stream dimensions obtained in this ap- 
proximation differ from (2) and (3) by a factor of 

\J zpi, however, analysis of the density distribution 
across the stream shows that the stream profile con- 
sists of a dense core whose characteristic dimensions 
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are close to expressions (2) and (3) and a rarefied pe- 
ripheral region. 

The dimensions of the stream are of fundamental 
importance for determining the characteristics of the 
stream of matter flowing from the vicinity of L\. The 
values for the total mass flow from the system in var- 
ious approaches are essentially equal (within a factor 
of two), and can be determined from the relation 

M - poCoS, (4) 

where p and Co are the gas density and sound speed 
at the surface of the outflowing star and 

j £.2 £,2 

S = ivab = 7iy g y (fj,)g z (fi)-^ ~ 

is the area of the stream cross section. Since for a 
( c 2 \^ 

polytropic gas p — I j , the final expression 
for the dependence of the mass flow on c is: 

If we transform the "energetic" overflow of the Roche 
lobe to " geometrical" overflow using the formula 

c% = A$ = -^r-AiZ, 

is the effective radius of the donor star; AR = 
R\ — Rlu where Rl x is the radius of a sphere with 
volume equal to the volume of the Roche lobe), then 
for 7 = 5/3, we obtain for the mass flow the widely- 
used expression 

M ~ (AR) 3 . 

These expressions for the mass flow into the system 
were derived in the framework of simplified models. 
Their common and, obviously most serious, disadvan- 
tage is the lack of allowance for the system's common 
envelope, which appreciably influences the flow pat- 
tern along the Roche surface and, correspondingly, 
the mass exchange parameters in the system. This ef- 
fect has unfortunately previously been neglected, due 
to the fact that a correct study of this phenomenon is 
possible only using three-dimensional numerical gas 
dynamics models. 



3. Numerical Simulation Results 

We base our analysis of the flow pattern near the 
donor star filling its Roche lobe on the results de- 
scribed in Paper I. The formulation of the problem 
is laid out in detail in Paper I, so that we repeat 
here only the basic features of the model, referring 
the reader to Paper I for details. 

(1) We considered a low-mass X-ray binary simi- 
lar to X1822-371 in which the mass of the outflowing 
component Mi = 0.28 M Q , the temperature of the 
gas on the surface of this component T — 10 4 K, the 
mass of the secondary (a compact object with radius 
O.O5i? ) - M 2 = 1.4M Q , the system orbital period 
Porb = l d -78 , and the distance between the compo- 
nent centers A = 7.35i?©. 

(2) A three-dimensional system of gas dynamics 
equations closed by the ideal gas equation of state 
was used to describe the flow. 

(3) To take into account radiative losses in the sys- 
tem, we took the adiabatic index to be close to unity, 
namely 7 = 1.01, which is close to the isothermal 
value [20]. 

(4) We assumed that the outflowing star fills its 
Roche lobe and that the velocity of the gas at the 
surface is directed normal to the surface. We took 
the gas velocity to be equal to the local sound speed 
u = c = 9km/s. The density at the surface of the 
outflowing component p was taken to be p a . Note that 
the boundary value of the density does not affect the 
solution, due to the scaling of the system of equations 
in p and P. We used an arbitrary value for po in the 
calculations; when a specific system with a known 
rate of mass loss is considered, the real densities in 
the system can be determined simply by adjusting the 
calculated density values using the scale indicated by 
the ratio of the real and simulated densities at the 
surface of the outflowing component. 

(5) We adopted the conditions of free mass outflow 
both at the accretor and at the outer boundary of the 
calculation region. 

(6) The calculated region was the parallelepiped 
(-A..2A) x (-A..A) x (0..A); by virtue of the prob- 
lem's symmetry relative to the equatorial plane, the 
calculations were performed only for the upper half- 
space. 

(7) We used a high-order Total Variation Dimin- 
ishing scheme to solve the system of equations on a 
non-uniform (finer along the line connecting the com- 
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ponent centers) difference grid with 78 x 69 x 35 nodes. 

The system of equations was solved starting from 
arbitrarily chosen initial conditions up to the estab- 
lishment of a stationary flow regime. To ensure that 
this flow regime was established, we continued the cal- 
culations over a rather long time interval (exceeding 
ten orbital periods) after the initiation of the station- 
ary regime. 

The overall flow pattern in the equatorial plane of 
the system considered is presented in Fig. 1, which 
depicts density and velocity field contours in the area 
from -5 to 13i?0 along the X axis and from -6 to 6Rq 
along the Y axis. Figure 1 also shows four stream lines 
(marked with the letters A, B, C, D), which illustrate 
the directions of flows in the system. Analysis of the 
results in Fig. 1 shows that a significant fraction of 
the common envelope gas (stream lines A and B) ap- 
proaches the surface of the outflowing star (Roche 
lobe) in the course of its motion, and is accreted by 
this surface, preventing the formation of flows along 
the stellar surface. Some of the envelope material 
(stream lines C and D) approaches to the donor star 
and, stripping gas off its surface, participates in the 
formation of the stream in the vicinity of L\. 

Details of the flow pattern in the vicinity of the in- 
ner Lagrange point are shown in Fig. 2, which shows 
the same flow parameters as in Fig. 1 in a small re- 
gion of the equatorial plane from -1 to 5i? along the 
X axis and from -2 to 2i? Q along the Y axis. We can 
see from Fig. 2 that the common envelope gas sub- 
stantially changes the flow structure in the vicinity of 
Li, in particular, it strips matter off the surface of 
the donor star. The asymmetry of the action of the 
envelope gas on the outflowing stream is also visible 
in Fig. 2. Gas flowing in from above, along the direc- 
tion of the orbital motion, is accreted by the donor 
star, and only in a small region near L\ is this mat- 
ter taken up from the surface and transported to the 
stream. The envelope gas, approaching to L\ from 
below, strips matter off a significant fraction of the 
donor star surface. 

This effect of " stripping" matter off the surface of 
the donor star substantially changes the generally ac- 
cepted views of the formation of the matter stream 
flowing into the system and of the corresponding 
mass exchange parameters. According to the stan- 
dard models, the structure of the atmosphere near the 
surface of the donor star is determined by the equa- 
tions of hydrostatic equilibrium (see, e.g., [1, 18]). For 
the adopted system parameters and gas temperature 
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Fig. 1. — Density isolines and velocity vectors in the 
equatorial plane of the system. Roche equipoten- 
tials are shown with dashed lines. Four stream lines 
marked with the letters A, B, C, D are also shown, 
illustrating directions of the flow of the common enve- 
lope material. The location of the accretor is marked 
with the dark circle. The vector in the upper right 
corner corresponds to a velocity of 800 km/s. 




Fig. 2. — Density isolines and velocity vectors in the 
equatorial plane of the system in the vicinity of the in- 
ner Lagrange point. Roche cquipotcntials are shown 
with dashed lines. The vector in the upper right cor- 
ner corresponds to a velocity of 300 km/s. 
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(sound speed) at the donor star surface, the gas en- 
ergy is not large enough for the gas to escape directly 
from the stellar surface, so that the mass flow asso- 
ciated with thermal escape will be negligible small 
compared to the flow through the vicinity of L\. Gas 
located in the near-surface layer (with characteristic 
height of the order of the atmospheric scale height) 
can flow along the stellar surface, but, in this case, 
also, the total mass flow into the system does not 
change significantly [1] . 

The situation fundamentally changes only when 
the system's common envelope is taken into account, 
because, in this case, gas from the surface layer can 
be "stripped" from the star and carried into the sys- 
tem. Our calculations have shown that this is the flow 
regime that is realized in the system considered here. 
The momentum of the rarefied common envelope gas 
along the surface of the star is large enough to pick 
up matter from the surface layer over a significant 
fraction of the stellar surface. Further, this material 
and the gas flowing out from the vicinity of L\ form 
the stream of matter flowing into the system, and also 
determine the total mass flow. 

The effect of the envelope gas on the flow structure 
near L\ is especially visible in Fig. 3, which shows 
a three-dimensional projection of stream lines orig- 
inating on the stellar surface and entering the sys- 
tem with the gas stream. The area containing the 
points through which the stream lines enter the sys- 
tem represents that part of the surface providing mass 
exchange in the system. We can clearly see the asym- 
metry of this area in Fig. 4, which also shows stream 
lines leaving the stellar surface (also in 3-D projec- 
tion), but for an observer located on the line join- 
ing the component centers. The size of this area 
substantially exceeds the size of the stream in stan- 
dard models (shown by the ellipses near L\ in Figs. 3 
and 4). Analysis of these results indicates that the 
stream forms in a region is more than order of mag- 
nitude larger than predicted by theoretical estimates 
obtained without allowance for the common envelope 
and with the same values for the gas parameters at 
the donor star surface. The matter flow injected into 
the system can now be represented as the sum of two 
terms, and, accordingly, we obtain instead of (4) 

M = M Ll + Msnvel, (5) 

where the first term, derived from standard estimates, 
is an order of magnitude smaller than the second 




Fig. 3. — Matter stream lines originating from the 
donor star surface and forming the gas stream. The 
stream cross-section obtained in standard models is 
shown by the solid ellipse near L\. 



v M 

Fig. 4. — View of the matter stream lines in Fig. 3 
for an observer on the line connecting the centers of 
the two stars. The stream cross-section obtained in 
standard models is shown by the ellipse near L\. The 
projection of the outflowing star on the plane depicted 
is shown by the outer ellipse. 
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term, determined by the common envelope gas. 
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Fig. 5a. — Density isolines and velocity vectors in 
the stream section formed by the YZ plane passing 
through the inner Lagrange point. The stream cross- 
section obtained in standard models is shown by the 
bold ellipse. The vector in the upper right corner 
corresponds to a velocity of 200 km/s. The maxi- 
mum density (at the stream center) is l-2po, and the 
minimum density (outer contour) is O.Olpo- 

We will now consider the distributions of gas pa- 
rameters across the stream cross section. Figure 5 
shows density and velocity vector contours in the 
stream cross section formed by the YZ plane pass- 
ing through L\. The maximum density (in the cen- 
ter) exceeds the density at the stellar surface po by a 
factor of 1.2; the minimum density, depicted by the 
outer contour, corresponds to O.Olpo- We can see from 
Fig. 5 that the matter stripped off the stellar sur- 
face by the common envelope gas collects in a small 
zone centered on the equatorial plane. The center of 
this zone, which is determined by the common en- 
velope gas flows considered above, does not coincide 
with L\ (the point (0,0) in the figure); moreover, it 
lies beyond the matter stream inferred from standard 
models (shown by the ellipse near L\ in Fig. 5). The 
two-humped nature of the density distribution can 
be seen particularly well in Fig. 6, which shows the 
density distribution (normalized to its value at the 
stellar surface) along the Y axis for the section of the 
stream under consideration. The stream boundaries 




Y/R 

Fig. 5b. — Density distribution (normalized to the 
value at the stellar surface) along the Y axis for the 
stream cross section presented in Fig. 5a. The stream 
boundaries obtained in standard models are marked 
with vertical lines. 

obtained in standard models are shown in the fig- 
ure by the vertical lines. The effects described above 
can be interpreted as the presence of two well-defined 
streams in the system, one associated with the gas 
flowing from L\, and the other forming under the ac- 
tion of the common envelope gas. Further analysis of 
the calculation results indicates that the distributions 
of gas parameters across the stream change with dis- 
tance from the inner Lagrange point, and the stream 
becomes a single flow at a distance of the order of 
0.1i? Q (0.013A) from L\. This picture is supported 
by the results in Fig. 6, which shows density contours 
in the Y Z cross section of the stream O.li?0 from 
the inner Lagrangian point. The maximum density 
in Fig. 6 is 1.4p , and the minimum (outer contour) 
is O.Olpo- 

4. Conclusion 

Analysis of the three-dimensional numerical sim- 
ulation results presented here indicates that taking 
into account the impact of the common envelope in 
semidetached binary systems radically changes the 
mass exchange parameters obtained. For the low- 
mass X-ray binary X1822-371 that we have consid- 
ered, the overall mass inflow to the system increased 
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Fig. 6. — Density isolincs and velocity vectors in the 
stream section formed by the Y Z plane passing O.li?© 
from the inner Lagrange point. The vector in the 
upper right corner corresponds to a velocity of 200 
km/s. The maximum density (at the stream center) 
is 1.4/9o, and the minimum density (outer contour) is 
O.Olpo. 



by a factor of about 50 compared to the estimates 
of standard models for the same gas parameters val- 
ues at the surface of the donor star. In addition, the 
common envelope gas changed the flow pattern near 
the surface of the outflowing component, which ulti- 
mately affected the overall structure of the gas flows 
in the system, and, thus, the expected observational 
manifestations of these flows. 

Unfortunately, the quantitative estimates we have 
obtained are valid only for the specific object consid- 
ered, and do not allow us to draw more general conclu- 
sions about changes in the mass exchange parameters 
expected for other semidetached systems. This is due 
to the fact that the mass flow increase depends not 
only on the parameters of the binary system, but also 
on the properties of the common envelope gas, which 
can only be determined using three-dimensional nu- 
merical simulations. Nevertheless, given the results 
presented here, it is clear that an adequate descrip- 
tion of the mass flows in any semidetached system 
with a common envelope is possible only in models 
that take into consideration the influence of this en- 
velope. 
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